NPS  ARCHIVE 
1962.06 
BURNS,  R. 


AN  INVESTIGATION  OF  THE  COUPLING  BETWEEN  LONGITUDINAL 

AND  FLEXURAL  MODES  OF  VIBRATION  CAUSED  BY  AN  OFFSET 
OF  THE  CENTER  OF  GRAVITY  FROM  THE  CENTER  OF  STIFFNESS 


by 


RICHARD  F.  BURNS,  LIEUTENANT,  U.S.  NAVY 
B.S.,  UNITED  STATES  NAVAL  ACADEMY 

(195*0 


SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE  RETIREMENTS 
FOR  THE  DEGREE  OF 

NAVAL  ENGINEER 

AND  THE  DEGREE  OF 

MASTER  OF  SCIENCE  IN  NAVAL  ARCHITECTURE 

AND  MARINE  ENGINEERING 

at  tlie 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

June  1962 


Thesis 
B88416 


Professor  Frank  M.  Lewis 
Thesis  Supervisor 


OUOLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  3CHOOI 
MONTPJ*EY  CA  93943-51  Dl 


Of  xl 


AN  INVESTIGATION  OF  THE  COUPLING  BETWEEN  LONGITUDINAL 

AND  FLEXURAL  MODES  OF  VIBRATION  CAUSED  BY  AN  OFFSET 
OF  THE  CENTER  OF  GRAVITY  FROM  THE  CENTER  OF  STIFFNESS 


fey 


RICHARD  F.  BURNS,  LIEUTENANT,  U.S.  NAVY 
B.S.,  UNITED  STATES  NAVAL  ACADEMY 

(195*0 


SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE  RETIREMENTS 
FOR  THE  DEGREE  OF 

NAVAL  ENGINEER 
AND  THE  DEGREE  OF 
MASTER  OF  SCIENCE  IN  NAVAL  ARCHITECTURE 
AND  MARINE  ENGINEERING 
at  the 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
June  1962 


Signature  of  Author  . 


Department  of  Naval  Architecture  and 
Marine  Engineering,  May  18,  1962 


Certified  by 


Frank  M.  Lewis,  Thesis  Supervisor 


Accepted  by  . 


Chairman,  Departmental  Committee 
on  Graduate  Students 


AN  INVESTIGATION  OF  THE  COUPLING  BETWEEN  LONGITUDINAL 

AND  FLEXURAL  MODES  OF  VIBRATION  CAUSED  BY  AN  OFFSET 
OF  THE  CENTER  OF  GRAVITY  FROM  THE  CENTER  OF  STIFFNESS 

by 
Richard  F.  Burns,  Lt.  USN 

Submitted  to  the  Department  of  Naval  Architecture  and  Marine  Engineering  on 
18  May  1962  in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Naval  Engineer  and  the  degree  of  Master  of  Science  in  Naval  Architecture  and 
Marine  Engineering. 

ABSTRACT 

The  purpose  of  the  investigation  is  to  find  how  the  natural 
frequency  of  transverse  vibration  is  affected  when  the  center  of  gravity 
and  center  of  stiffness  are  not  coincident  and  to  determine  if  it  is 
important  in  calculating  the  natural  frequency  of  the  fundamental  mode 
of  vibration  of  ships. 

Two  methods,  one  an  approximate  analytic  method,  and  the  other 
a  lumped  mass  method,  are  used  to  determine  the  strength  of  the  coupling 
of  these  modes  of  vibration.  The  approximate  analytic  method  involves 
the  solution  of  differential  equations  from  an  assumed  mode  shape  and 
then  the  use  of  Stodola's  equation  to  compute  the  change  in  natural 
frequency .  The  lumped  mass  method  is  a  numerical  procedure  requiring  a 
considerable  amount  of  computation.  Both  methods  are  in  good  agreement 
when  the  ratio  of  natural  transverse  to  longitudinal  frequency  is  less 
than  one-half. 

The  approximate  analytic  method  shows  that  the  change  in 
transverse  frequency  of  vibration  is  a  function  of  the  product  of  two 
factors.  One  factor  is  directly  proportional  to  the  area  times  the  distance 
of  the  center  of  gravity  from  the  center  of  stiffness  squared  and  inversely 
proportional  to  the  area  moment  of  inertia.  The  other  factor  is  a  com- 
plicated expression,  but  which  is  dependent  only  on  the  ratio  of  trans- 
verse to  longitudinal  frequency  as  computed  by  simpler  theory  in  beams, 
neglecting  the  offset  of  the  center  of  gravity  from  the  center  of  stiffness. 
These  expressions  have  been  evaluated  and  plotted  against  frequency  ratio. 
A  simple  method  of  computing  the  change  in  transverse  frequency  of  vibration 
is  provided. 

When  the  important  parameters  for  calculating  the  change  in 
frequency  of  vibration  are  calculated  for  a  merchant  ship  and  submarine, 
it  is  seen  that  the  coupling  of  the  fundamental  modes  of  vibration  in 
these  structures  is  insignificant. 

In  the  higher  modes  of  vibration,  the  ratio  of  transverse  to 
longitudinal  frequencies  approaches  unity.  It  is  recommended  that  the 
coupling  in  these  higher  modes  be  investigated. 
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NOTATION 

A      area  of  a  cross  section  which  will  support  longitudinal  stress 

c      arbitrary  constant 

E      modulus  of  elasticity 

f      ratio  of  transverse  to  longitudinal  frequency  computed  disregarding 
the  offset  of  the  center  of  gravity  from  the  center  of  stiffness 

F  force 

g  gravity 

I  area  moment  of  inertia  of  a  cross  section 

K  spring  constant 

L  length  of  a  beam 

I  distance  between  masses  in  the  lumped  mass  method 

M  moment 

m  mass 

5  stress 

u      longitudinal  displacement 

V      shear 

w      frequency 

w      longitudinal  frequency  as  computed  disregarding  the  offset  of  the  center 
of  gravity  from  the  center  of  stiffness. 

w  .     transverse  frequency  as  computed  disregarding  the  offset  of  the 
center  of  gravity  from  the  center  of  stiffness. 

w,  transverse  frequency  with  coupled  longitudinal  motion 

y  vertical  displacement 

y  distance  of  the  center  of  gravity  from  the  center  of  stiffness 

6  slope 

XL     a  number  determined  by  the  transverse  mode  of  vibration.  It  equals 
4.73  for  the  fundamental  mode. 

(ctA)    weight  per  unit  length 
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HSTROIXJCTION 

The  calculation  of  the  natural  frequencies  of  the  longitudinal  and 
transverse  modes  of  vibration  of  any  structure  have  long  been  considered 
independent.   '  In  a  uniform  bar,  of  course,  it  is  true,  but  for  example, 
in  a  ship's  hull  where  the  center  of  mass  and  center  of  stiffness  are  not 
coincident,  it  is  not  necessarily  true.  A  bending  of  the  ship's  structure 
causes  some  motion  in  the  longitudinal  direction  of  the  center  of  gravity. 
This,  in  turn,  will  cause  some  change  in  the  natural  frequency  of  transverse 
vibration.  It  is  the  purpose  of  this  thesis  to  investigate  the  strength  of 
the  coupling  of  longitudinal  and  transverse  modes  of  vibration  when  the 
centers  of  gravity  and  stiffness  are  displaced  and  determine  if  it  has  any 
significance  in  the  calculation  of  natural  frequencies  of  ships. 
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PROCEDUEE 

Two  distinct  procedures  are  used  to  find  the  coupling  of  longitudinal  and 
flexural  modes  of  vibration;  one  method  is  an  approximate  analytic  method 
and  the  other  is  a  lumped  mass  method.  The  steps  to  each  method  are  outlined 
in  subsections  A  and  B  below;  the  details  are  in  the  appendices.  Through 
a  numerical  example  the  results  of  the  two  procedures  are  then  compared 
for  agreement.  The  entire  analysis  is  then  interpreted  in  the  light  of 
typical  values  of  significant  parameters  for  a  submarine  and  a  merchant  ship. 

The  investigation  is  restricted  to  the  fundamental  mode  of 
vibration. 

In  neither  the  approximate  analytic  method  nor  the  lumped  mass 
method  is  the  deflection  of  the  system  due  to  shear  and  rotary  inertia 
considered.  The  transverse  motion  is  due  to  bending  only. 
A.  Approximate  Analytic  Method. 

The  general  procedure  in  this  method  is  to  find  a  correction  to  an 
assumed  mode  shape  from  the  differential  equations  of  a  system  in  which 
bending  causes  longitudinal  motion  and  longitudinal  forces  cause  bending. 
The  change  in  mode  shape  is  then  used  as  a  measure  of  the  change  in  natural 
frequency  of  transverse  vibration.  The  steps  in  the  procedure  which  are 
carried  out  in  Appendix  I  are; 

1.  Determine  the  differential  equations  of  motion  of  a  system  in 
which  the  center  of  gravity  and  center  of  stiffness  are  not 
coincident.  The  first  equation  is  a  longitudinal  force  balance 
equation  and  the  second  one  is  a  bending  moment  equation. 

.  Mi  i-  (u  +  |z  y)  +  ae  A  .  0  (1) 

8   dt2     8x      ax2 
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(|£>^  +  EI^  +  EA^u  -  =  0  (2) 

6   dtc       dx*     dx^ 

(aA)  =  weight  per  unit  length 

A    =  area  of  cross  section 

g    =  gravity 

u    =  longitudinal  displacement 

y    =  vertical  displacement 

E    =  Young's  modulus 

y    =  distance  of  center  of  gravity  from  the  center  of 
stiffness 

I    =  area  moment  of  inertia  of  a  cross  section 

2)  Assume  a  transverse  mode  shape  of  a  free-free  beam  in  space. 

'  -  c[cos  x  *  +  Cn^f  >  Cosh  x  x  ]  <3) 

c  =  arbitrary  constant 

XL  =  a  number  determined  by  the  mode  of  vibration.  It  is 
equal  to  4.73  for  the  fundamental  mode  of  vibration. 

L  =  length  of  the  structure. 

3)  From  the  differential  equation  for  longitudinal  force  balance, 
find  the  longitudinal  motion  of  the  center  of  stiffness  pro- 
duced by  the  mode  shape. 

k)     This  longitudinal  motion  causes  longitudinal  forces  which 

produce  moments  about  the  center  of  gravity  and  hence  additional 
bending  of  the  system.  Determine  the  additional  deflection 
from  these  forces  using  a  bending  moment  equation. 

5)  With  the  results  of  Step  k   recycle  Steps  3  and  k.     Hence,  a 
first  and  a  second  approximation  to  the  additional  deflection 
can  be  obtained. 


6)  Use  a  method  similar  to  that  of  Stodola  to  find  the  trans- 
verse frequency  with  coupled  longitudinal  motion.  The  ratio 
of  the  square  of  the  transverse  frequency  with  longitudinal 
motion  to  that  without  is  proportional  to  the  ratio  of  total 
deflection  without  longitudinal  motion  to  the  total  de- 
flection of  the  beam  with  longitudinal  motion. 

B.  Lumped  Mass  Method. 

The  lumped  mass  procedure  is  the  Prohl-Myklestad  sequence  method 
for  calculating  the  natural  frequency  of  beams  of  variable  cross  section. 
This  calculation  is  modified  to  include  the  effect  of  the  displacement  of 
the  center  of  gravity  from  the  center  of  stiffness. 

The  calculation  is  based  on  the  fundamental  equations  of  bending 
and  strain  in  beams: 

lL  =  il  and 

EI   .  2 
dx 

t 
The  beam  is  divided  into  sections.  At  each  section  an  inertia  loading  is 
added? 

in  shear  equal  to 

2 

m  w  y 

in  bending  equal  to 
m  w  y  (u  -  y  e) 

and  in  longitudinal  force  equal  to 
m  w  (u  -  y  9) 

m  =  mass 

w  ■  frequency 

9  =  slope 

I   =  length 
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The  procedure  is  to  calculate  the  longitudinal  force  and  displacement,  trans- 
verse deflection,  slope,  bending  moment  and  shear  at  each  section  along  the 
beam  in  proceeding  from  the  end  of  the  beam  to  the  mid-point  <>  Various  fre- 
quencies are  tried  until  the  boundary  conditions  are  satisfied.  For  the 
fundamental  mode,  the  slope,  shear,  and  longitudinal  motion  at  the  mid-point 
are  zero. 
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EESULTS 


A.  Approximate  Analytic  Method 

The  equation  for  determining  the  frequency  of  transverse  vibration 


by  the  approximate  analytic  method  is: 

Y 
2  _         2  o 

Wt     =     ot       Y0  +  Y1+Y2 


W 


where : 


w 


ot 


Egl 


(5) 


YQ  =  c[2  Cos      (XL/2)  -  1  - 


A     -2 

Yi  =  T   y    c  Cos  *L' 


p  k       Mi 

»2  D.     X      (Cos|Dj~  -  1) 


(6) 


(Cos  XL/2  -  1)  1     (Cosh  XL/2 

(1  -  D1  X2)  Cos  XL/2     (1  +  D.^2)  Cosh  ; 


(7) 


2     k         ,    2     k 
.  2   .  r  =2  D  *  X*    k  D  *  X* 

.  (T)  y  c  cos  W*  [(—^  -  (i  ■  -i  xV)- 
c°°^  -1     d,2^     L/2 

( — Its—)  -  0     ,  l)  (*£•> 


1  -  D 


i  X    ^'   (Cos^f) 


[Cos  XL/2  -"!)      j. 
(1-^X2)2    C0S  XL/2       U^D^2)2' 


(Cosh  XL/2  -1)1 
Cosh  XL/2   J 


(8) 


w.  =transverse  frequency  with  coupled  longitudinal  motion 

w  .-transverse  frequency  as  calculated  disregarding  the  offset 
°  of  the  center  of  gravity  from  the  center  of  stiffness 


Y  =  total  deflection  of  the  assumed  mode  shape 
o 
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Y  =  first  approximation  to  the  additional  deflection 
Y«  -  second  approximation  to  the  additional  deflection 

When  substitutions  shown  in  Appendix  A  are  made,  all  the  terms  in  brackets 
in  equations  (7)  and  (8)  are,  to  an  approximation,  dependent  only  on  fre- 
quency ratio,  f . 

wo£ 
where 

w   =  longitudinal  frequency  as  calculated  disregarding  the 

offset  of  the  center  of  gravity  from  the  center  of  stiffness.  Specifically, 

t   ot 


Similarly, 

^   2 

Another  important  result  derived  in  Appendix  A  is  that  the  frequency  ratio, 

f,  is  independent  of  mass  per  unit  length. 

These  two  important  results  suggest  a  rearrangement  of  the 

frequency  equation  h,   to  permit  easy  calculation  of  the  change  in  frequency 

for  any  system. 
2 

£-)  -      x 


w  '       Y    Y 
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This  equation  can  "be  put  in  a  still  more  usable  form. 
2 


(#-) 


i  *  €■   s(f)  +(^)2  Mf ) 


(9) 


I   ««-'  ■*  i 

where  g(f )  and  h(f ),  functions  which  are  only  dependent  on  frequency  ratio, 
can  be  obtained  from  Figures  I  and  II. 

The  function  for  the  first  approximation  to  beam  deflection  is: 


Cos 


XL 


g(f)  = 


h 

I  (^)  (Cos  f?  -  1) 

■rm^fn'  %  2         '         (Cos  XL/2 

%  '  {i  f)4)Cos  (fi}  (1  "(^)2)  Cos  ^     (1  +  (^)2) 


,  (Cos  XL/2  -  1) 
(1  +   (— ' 

2  V  VfTT' 


\z  Cos  XL/2  -  1  -  g'VffiJ 
L  '  Cosh  XL/2J 


The  function  for  the  second  approximation  to  beam  deflection  is: 


k  k 

2(— )     4  (—) 


,XL, 


Cos  f  5  -  1    (~) 
h(f )  =  Cos  XL/2  |(— ~5  -      XL  4  2><     tr    '  >"  (J^S 

L  1  -  (~r       (1  -  (~)  f  Cos  f|        1  -  (^)2 


) 


» .  (fZ)  sin  f  2  +  (cos  xl/2  -  i; 

2   (Cos  f?)2   (1  -  (^)2)2  Cos  XL/2 


'fir 


•* 


- |  i  [  2  Cos  XL/2 

(1  +  (^)  )  Cosh  XL/2  J  '  L 


Cos  XL/2 

"xl72 
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These  expressions  do  not  have  to  be  evaluated  in  finding  the  change  in 

Ay2 
natural  frequency.  All  that  has  to  be  evaluated  is  w  .  ,  w  ,  and  h£-  with 

g(f )  and  h(f )  from  Figures  I  and  II  to  use  the  approximate  analytic  method. 

Above  a  frequency  ratio  of  one-half  the  approximate  analytic  method 

gives  poor  results.  However,  it  is  still  desirable  to  estimate  frequency 

changes  at  higher  frequency  ratios.  An  examination  of  the  relative  order 

of  magnitude  of  these  two  functions  suggests  an  alternate  method  of  determining 

frequency  changes  from  data  of  only  one  numerical  example  of  the  lumped  mass 
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FIGURE  I 


FUNCTION  FOR  THE  FIRST  APPROXIMATION 
TO  THE  CHANGE  IN  BEAM  DEFLECTION 


g(f )  versus  f 
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FIGUEE  II 


FUNCTION  FOR  TEE  SECOND  APPROXIMATION 
TO  THE  CHANGE  DT  BEAM  DEFLECTION 
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method  shown  in  Figure  IV.  This  alternate  method  is  based  on  the  fact  that 
h(f )  «  [g(f )] 

To  use  Figure  IV  for  other  beams  use  the  following  procedure: 

wt 

1.  Determine  —  from  this  figure  for  the  desired  frequency  ratio, 

ot 

2.  This  number  can  be  used  to  solve  Equation  (10)  for  x.  The 
solution  is  shown  as  Equation  (ll). 

3.  Use  this  value  of  x  in  a  new  frequency  equation  with 


-2 

I 


for  the  desired  application.     See  Equation  (12). 
2 

<T  }    = " 2-2  (10) 

ot  1  +  ax  +  a      x 

o  o 


2 
w 


O  *  X 

w-  2 

tr->  " 4 — 2-2  <12> 

ot  1  +  a.,   x     +  a,     x 


where 

wt 

— -  =  frequency  change  obtained  from  Figure  IV 

wot 

-2 

A?2 
a,  =  — i-  for  the  desired  application 

—  =  frequency  change  for  the  desired  application. 
Wot 

B.  Results  of  the  Lumped  Mass  Method. 

The  equations  for  determining  the  conditions  at  each  mass  along 
the  beam  are  shown  below.  The  zero  subscript  refers  to  the  conditions  at  the 
previous  section;  the  1  subscript  refers  to  the  mass  at  the  point  being 
evaluated.  The  actual  numerical  value  of  each  of  these  equations  becomes  the 
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LUMPED  MASS  SYSTEMS 
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subscript  zero  term  when  the  next  point  along  the  beam  is  evaluated. 
Ml  =  Mo  +  V  +  ^  *1  w2  ^ul  "  y  91^ 

Vl  "  Vo  +  "l  w2  yl 

Mi       V  Z2 
91  "   9o  +  EI     +  2EI 

2  3 

m  r     vr 

yi =  yo +  V  +  ar +  mr~ 

F 
o 

1         o      K 


F.,    ■  F 


'0  +  *!  w    (^  -  y  e2) 


where 


M  =  bending  moment 

V  =  shear  force 

F  =  longitudinal  force 

K  =  spring  constant  between  masses 

i  =  distance  between  masses 
One  advantage  of  the  lumped  mass  method  is  that  it  is  useful  in  deriving  a 
frequency  equation.  For  the  three  mass  systems  shown  in  Figure  III,  the 
frequency  equation  isc. 

n   .  9      m     2  ,1   £3  2  I  y2  Aw  2s         2  k  ,tk  y2  .  1  X3  x 

C.  Results  of  a  Numerical  Example. 

In  order  to  check  the  validity  of  the  results,  a  numerical  example 
was  done  using  both  methods.  The  change  in  frequency  is  shown  plotted  in 
Figure  IV  as  a  function  of  frequency  ratio,  f .  The  plot  is  for  a  system  of 
these  properties: 
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A  .  1  in2 

I  =  1/12  in 

y  =  l/k  in 

a  -  .284  lbs/in3 

*£    -75 

In  the  lumped  mass  method  the  system  was  divided  into  eleven  masses  (Figure  III) 
For  system  of  fewer  masses  the  agreement  between  methods  was  poor.  For  example, 
at  a  frequency  ratio  of  one-fourth: 


wot 


=  .974   analytic  method 


t 
— -  =  .961   eleven  mass  system 


Ot 


Wot 


Wot 


=  .950   nine  mass  system 


=  .910   seven  mass  system 


For  all  numerical  work,  w   and  w   were  computed  from: 


w   =  ^ 
Wot    L2 


22A   /gEI 

L2    V^ 


ot       L   V  TaAj 

D.  Numerical  values  for  Marine  Application. 

The  approximate  analytic  method  indicates  that  there  are  two  sig- 
nificant parameters  in  calculating  the  change  in  transverse  frequency,  namely, 

Ay2 
frequency  ratio,  f,  and  the  geometric  factor  -*r-  .  The  geometric  factor 

has  been  evaluated  for  two  ships: 


-16- 

2 
merchant  ship  (3)    =SL  «  .087 

2 
submarine  (4)       -£-  «*  .017 

These  numbers,  of  course,  are  dependent  on  the  condition  of  loading  and, 
furthermore,  they  are  valid  only  for  the  midship  section.  The  geometric  factor 
for  the  merchant  represents  an  upper  limit  of  this  number  for  the  "Gopher 
Mariner".  Approximate  frequency  ratios  for  these  ship  types  ares 

submarine  (k)         f  ^  .15 

merchant  ship    f  m  .25 
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DISCUSSION  OF  RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS 

A.  Accuracy  of  Methods 

The  approximate  analytic  method  is  inaccurate  above  the  frequency 
ratio  of  one-half.  The  lumped  mass  method,  on  the  other  hand,  can  be  used 
over  the  entire  range  of  frequency  ratio. 

Reference  (5)  states  that  to  ensure  better  than  one  percent 
accuracy  in  computing  frequencies  by  the  lumped  mass  method,  thirteen 
masses  should  be  used  per  wave  length.  If  more  masses  had  been  used  in  the 
numerical  example,  the  agreement  between  methods  at  the  low  frequency  ratio 
would  be  better  than  it  already  is.  For  eleven  masses  there  is  only  a  1.3$ 
difference . 

B.  Significance  in  Ship  Vibrations. 

When  the  constants  for  the  submarine  and  merchant  ship  stated  in 
Part  III  above  are  used  in  Equation  (9)  in  conjunction  with  Figures  I  and  II, 
the  change  in  transverse  frequency  is  truly  insignificant  for  the  fundamental  mode 
For  higher  modes  of  vibration,  the  frequency  ratio  is  larger  and  the  change 
in  frequency  might  not  be  so  insignificant. 

C .  General 

Two  general  methods  for  calculating  frequency  change  due  to  an 
offset  of  the  center  of  gravity  from  the  center  of  stiffness  are  given  . 
One  method  is  approximate;  the  other  is  exact.  Also,  a  method  of  estimating 
frequency  change  by  a  simple  calculation  from  one  lumped  mass  numerical 
example  is  explained.  The  coupling  of  the  transverse  and  longitudinal  modes 
of  vibration  in  any  structure  can  be  estimated  or  obtained  exactly  by  the 
methods  explained  in  the  text  of  this  thesis. 
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All  the  work  is  restricted  to  the  fundamental  mode  of  vibration, 
The  lumped  mass  method,  of  course,  could  be  used  for  higher  modes .  It  is 
recommended  that  the  coupling  due  to  the  offset  of  the  center  of  gravity 
from  the  center  of  stiffness  for  higher  modes  he  investigated. 
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APPENDIX 
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APPENDIX  A 
DETAILS  OF  THE  APPROXIMATE  ANALYTIC  METHOD 

The  approximate  analytic  method,  outlined  below,  is  divided  into 
steps,  the  numbering  of  which  corresponds  to  the  numbering  in  the  procedure 
Section  IIA  above.  The  description  which  appears  in  that  section  is  not 
repeated . 
Step  1:  Derivation  of  the  Equations  of  Motion 

With  the  sign  convention  shown  in  Figure  A-l  an  equilibrium 
force  balance  on  a  particle  of  the  system  yields  this  differential  equation 
of  motion. 

ot  dx 

where : 

(aA)  -  weight  per  unit  length 

g  =  gravity 

y  =  y  -  y  =  distance  of  center  of  stiffness  from  the  center  of 
gravity 

E  =  Young's  modulus 

u  =  displacement  in  the  x~  direction  of  the  center  of  stiffness 
from  the  equilibrium  position 

S  =  E#i 
dx 

ft2  2 

Letting       -  2-x  (      )     =    w*  (      ) 

at* 

the  equation  can  be  rewritten 

Ls*l  v2  (  u  +  ft  7)  ♦  EA^  «  0  (la) 

s  ^       dx* 

A  moment  balance  about  the  center  of  gravity  of  a  small  particle 

can  be  written  using  the  sign  convention  of  Figure  A- II. 
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figuhe  a-i 
longitudinal  force  balance 
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FIGUKS  A- II 

MOMENT  BAIANCE 
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Vdx  +  M-(M+dM)  +  (S+^|dx-S)Ay  =  0  (2) 

dV  JM)  il   dx 

g   dt2 


*L      d2y 
EI=  "^ 

du 

5  =  E  ^ 

dx 

A  =  area  of  the  cross  section 
Equation  (2)  can  be  rewritten; 

f    A    +   «i^  +  EA^y  =  0  (2a) 

6  dt*        dx        Sx"3 


Step  2:  Assumed  Mode  Shape  .  ,  i  J.  Koti  m 

According  to  the  classical  theory  oi'  vibration,  the  equation  of 
transverse  motion  for  a  uniform  beam  free  in  space  is: 

'  m  c[Cos  xx  +  (cosh\L02)  Cosh  H  <3) 

L  =  length  of  the  beam 

x  =  distance  from  mid- length 

c  =  arbitrary  constant 

For  the  fundamental  two=noded  mode  of  vibration: 


XL  -  4.73  =     n/  w2    gf  L  (3a) 

Step  3:     Longitudinal  Motion 

Equation  3  is  substituted  into  Equation  2  and  the  resulting 

differential  equation  solved  for  the  longitudinal  motion,  u.  The  arbitrary 

constants  are  evaluated  by  assuming  no  longitudinal  motion  at  mid-length 

and  no  stress  at  the  ends* 
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u     =  0 
du 


x  =  0 

x=L/2  g  =  0 

The  equation  for  longitudinal  motion  becomes: 


u  =  F.,   Sin  — -  +  F    Sin  Xx  +  P    Sinch  Xx 


where 


1         ^ 


f       .   B    XD     [l»i^]  [*>*  xl/2       ] 

1  X  2     Ll  -  Dx2  X4J  LCos  L/2^  J 


D    -«L 

w     (aA) 


D2  =  y  c     X 


Cos  XL/2 
3  '      "  Cosh  XL/2 


W 


(5) 

(6) 
(7) 


Step  hi     Additional  Deflection  Caused  by  Longitudinal  Forces 

In  this  step  the  objective  is  to  find  the  additional  vertical 
deflection  caused  by  the  longitudinal  motion  and  forces.  Again  referring 
to  Figure  A-2,  the  differential  equation  to  be  Satisfied  is: 
dM  .  ds 


ax  +  s  *ym0 


(8) 


where 


M_ 
EI 


dx 


r.  dU 

dx 


y'  =  additional  transverse  deflection 
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Equation  (8)  becomes 

&  -  -  £  7    4  0) 

dx3       x    dx 

These  boundary  conditions  are  applicable  in  the  solution  of  equation  (9); 

itf  -  0  x  -  L/2 

dx 

?£l  =  o  x  =  0 

dx 

u  =  0  x  =  0 

The  result  of  substituting  Equation  (k)   into  Equation  (9)  and 
solving  the  differential  equation  is  the  equation  for  the  additional  deflection 
of  the  beam  caused  by  the  coupled  longitudinal  motion.  The  equation  is; 

+ 
'1 

A  measure  of  the  total  additional  deflection  in  the  beam  is  obtained  by 

evaluating  Equation  (10 )  at  midlength  and  at  the  end. 

Yl=y'x=0  -  y'X  =  L/2  (11) 

,2  xk     (CoSi/^  -  1) 


y»  =  j    |Fi  \/Di  CoSi/5"  +  y~    CoS  Xx  "  X^  Cosh  Xxl  +  Constant  (10) 


Ay2,.   .  T  lo   r  2  D  1  X    — \/^T'     ,   Cos  XL/2  -  1 
=  ~--cCos  XL/2  —    p—r  — fTp---   ••  +  ■* — '    2 

Ll  -  D  d  X   Cos  yip-        (1  "  D!  x  )  Cos  XL/2 


^Cosh  XL/2 


.      .2   Cosh  XL/2 


(12) 


Step  5;  The  Second  Approximation  to  Total  Transverse  Deflection; 

When  Equation  (10 )  is  then  cycled  through  Steps  3  and  h  and 
evaluated  at  the  mid-length  and  end  point,  the  result  is  a  second  approximation 
to  additional  deflection  of  the  beam. 
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a2   -4  r"2Di2^        ^Di2^      Co#-1 

Y2  =  (§)     (y)%  Cos  XL/2  [(- i-g-Tf  \2,k){„         L/g       } 


1    -    D,       X  1    -    D,       X  COS     r=r- 

1  1  ^ 


2       If  r— 

(     "l  *   wL/2.  Sin  L/2\Tl  l  (Cos  XL/2  -  l) 

'    I3I5  ViT  (Cos  W     (1  -  Z  X2)2      Cos  XL>2 


+  i A(Cosh  XL/2  -  1)1         ,      ) 


Step  6;  Calculation  of  Natural  Frequency. 

With  an  assumed  frequency  and  deflection,  the  Stodola  procedure 
is  to  calculate  the  inertia  loading  and  resultant  deflection.  If  the 
assumed  frequency  equaled  the  natural  frequency,  the  resultant  deflection 
would  equal  the  assumed  deflection.  Mathematically ; 

wi2  *i  -  \  y0  (l4) 

w2  =  w2  -i  (15) 

1    o   y2 

y  =  assumed  deflection 
o 

y  =  resultant  deflection 

w  =  assumed  frequency 

w  =  natural  frequency  of  the  "beam 

When  the  Stodola  equation  is  applied  to  the  approximate  analytical 
procedure,  the  equation  for  determining  the  natural  frequency  of  a  system 
with  a  coupling  of  longitudinal  and  transverse  motion  becomes: 

w2  =  w  2 


t     ot   Y-*  +  Yc~+  y 
o    1    2 


Y  is  Equation  (3)  evaluated  at  midlength  and  at  the  end 

Y..  is  Equation  (12) 
J. 

Yp  is  Equation  (13) 
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w   =  natural  transverse  frequency  of  the  system  without  coupling 

ot 

of  longitudinal  and  transverse  motion, 
w   =  natural  frequency  of  the  system  with  coupling  of  longitudinal 
and  transverse  motion. 

Additional  information  can  be  obtained  by  some  algebraic 

manipulation.  From  Reference  (2),  the  equations  for  the  square  of  the 

natural  longitudinal  frequency  and  transverse  frequency  of  a  uniform  beam 

are  respectively: 

2 
lu     \2  -  (™\       EgA 


(Wot)  -TT"  loJT 


Therefore,       0 

f  -     2  -  ^r^  72   A  {lb} 

w  .      n     L 
oX 

Equation  (l6)  gives  the  very  important  result  that  the  frequency  ratio,  f, 

is  independent  of  mass  per  unit  length. 

p 
An  examination  of  the  expression  D.^  X  yields  another  important 

2 
result ,  Equations  (3a)  and  (5)  show  that  D,  X  is  equal  to: 

D  x* „  r^L_i  (itji3)2 

1  Lv?  (aA)J       L 


mL.     f£A     ±-      (4.73  r 

2   LToa7   t2  J  ^    f:u 

t 

t) 

If  the  change  in  transverse  frequency  is  small,  w^  w   , 
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Therefore, 

Similarly,  it  can  be  shown  that: 

4*         2 

These  manipulations  show  that  Equations  (12)  and  (13)  are  made  up  to  two 

Ay2 
factors;  one  is  -*r-  and  the  other  is  a  complicated  expression  which  is 

dependent  only  on  frequency  ratio. 
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AFPENDIX  B 
DETAILS  OF  THE  LUMPED  MASS  SYSTEM 

The  equations  in  the  beam  sequence  calculation  are  derived  from 
the  fundamental  equation  for  bending  in  a  beam. 

EI  =Jk   -  M 
dx 

In  the  case  of  a  beam  with  the  center  of  gravity  and  center  of  stiffness 

coincident,  this  equation  and  its  integrals,  using  the  sign  convention 

of  Figure  B-I,  became: 

,2 
EI  2JL   =  m  +  V  x 

d?    °    ° 

,  2 

EI  ?  =  M  x  +  V  ~  +  C- 
dx    o     o  2     1 

2      1 
Mx    V  xD 

Ely  =  -|—  +  -g—  +  c±x   +  c2 

where 

V  =  shear  at  the  boundary  or  beam  end 
o 

M  =  moment  at  the  boundary  or  beam  end 

When  these  quations  are  used  in  a  step-by- step  procedure,  the  equations 

for  moment,  shear,  slope,  and  deflection  when  evaluated  at  an  incremental 

distance  to  the  right  of  mass  1,  shown  in  Figure  B-I,  become: 

M  =  M     +  V  i 
o        o 

V  =  Vo   ♦  ^  W2  yx 

M     jt       v/ 

9  -  6o  +  ST"  +  2ET 

2  3 

Mi      vr 

y  -  yG  +  V +  A-  *  sir 
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FTGURE  B-I 

SIGN  CONVENTION  FOR  THE  LUMPED  MASS  SYSTEM 
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Center  of  Stiffness 
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where 

y  =  deflection  at  the  end 

9  =  slope  at  the  end 
o 

In  the  case  of  where  the  center  of  gravity  and  center  of  stiffness 

are  not  coincident,  these  equations  are  modified  and  two  additional  equations 

are  added. 

HL  =  Mo  +  VOyiD1w2(ury  e^  (l) 

Vl  "  Vo  +  "l  ^  yl  (2) 

Mi       V  I2 

en  ■  e    +  ~-  +  7&r-  (3) 

1         o       EI         2EI  VJ' 

Ml2      V  £3 

'i  -  *o +  V  +  sir +  sir  w 

F 

ui  =  u0-ir;  (5) 

Fi  =  Fo  +  "i  "~  ^ui  "  "~    i^  ^ 

where 

u  =  longitudinal  displacement  at  position  zero 

u  =  longitudinal  displacement  at  position  one 

AE 
K  =  spring  constant  =  — 

F  =  longitudinal  force  at  the  end  of  the  beam 
F,  =  longitudinal  force  at  position  1 

y  =  distance  of  the  center  of  gravity  from  the  center  of  stiffness 
In  the  case  of  a  beam  vibrating  free  in  space,  the  conditions  at  the  end 

of  the  beam  are  further  specified. 

2       - 

F  =m  w  (u  -y  9  ) 
o    o    x  o     o' 

M  =F  y=m  w  y(u  -y  8  ) 
o    o     o    ^yoo 

2 

V  =  m  w  y 
o    o    ^o 
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Equations  (l)  through  (6)  can  be  evaluated  at  other  positions  along  the  beam 
using  a  table  as  shown  in  Appendix  C»  To  start  the  calculation,  let 


yo 

=  1 

u 
o 

=  u 
o 

9o 

=  e 

o 

Then  a  table  is  completed  for  each  mass  along  the  system  using  the  results 
of  the  previous  table  to  start  the  next.  The  end  of  the  last  table,  which 
will  be  the  conditions  at  the  mid-point  of  the  beam,  will  yield  three 
equations; 

6  -  all  +  a21  uo  +  a31  6o 

V  =  a12  +  a22  uo  +  a32  9o 

u  "  a13  +  a23  U°  +  a33  8o 
If  the  proper  frequency  was  used  in  the  calculations,  these  equations  all 
equal  zero.  The  procedure  is  to  let  u  and  V  equal  zero,  solve  for  u  and 
9  and  substitute  into  the  e  -  equation.  When  6  equals  zero,  the  proper 
frequency  was  used  in  the  calculations. 
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APFENDIX  C 
SAMPLE  CALCULATIONS  OP  LUMPED  MASS  METHOD  AT  A  FREQUENCY  RATIO 
OF    .25  WITH  ELEVEN  MASSES 


A 

s  1  in. 

I 

=  1/12  in.4 

a 

=  .284  lbs. /in.3 

L 

=  8.24  in. 

f 

=  .25 

y 

=  1/4  ia- 

i 

=  .824  in. 

m 

=  .60666   x  10° 3 

( lbs.» sec  1 

.  /in.) 

2 

ot 

0 
=  3.69543  x  10 

( radians/ £ 

\2 
sec. ) 

-|  =  .03295  x  10"5  (l/lbs.-in) 

L--  *  .01357  x  10"5  (1/lbs.) 

^j  =  .00372  x  10°5  (In. /lbs.) 

i  =  .00274  x  10"5  (in. /lbs.) 


assume  w  1 

=  «975 

ot 

Therefore 

> 

2 

mw  = 

2.13120 

X 

105 

F       =  I.O656O  x  105  u     -    .26640  x  105   6 
00  c 

M      =   .26640  xlO^u     -   .06660  x  105  9 
000 

V   =  I.O65O  x  105  y 


o       '       vo 
The  five-page  tabular  calculation  to  reach  the  center  of  the  beam  is  on 
the  following  pages. 
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The  last  page  of  the  tabular  calculation  yields  three  equations  which  would 
all  equal  zero  if  the  right  frequency  was  assumed  at  the  beginning  of  the 
calculation. 

e  ■  1.27136  y  +  .23484  u  +  2.14740  e 

V  =  13.68628  y     +   .88648  u     +  23.551V7  9 

u  =  .00219  y  +  .92853  u  +  .01907  e 

let      yo  =  1 
u  =  0 

V  =  0 

Solve  for  u  and  9  and  substitute  into  the  9  equation. 
00 

uQ  =  .00958 

9  =-.581*48 

9  =  .02493 
When  the  assumed  frequency  is  .95  w  +>  &&&  the  five-page  tabular  calculation 
repeated, 

9  =  -  .01893 

The  result  is  plotted  in  Figure  C-l,and  thus  the  natural  transverse  frequency 

with  offset  center  of  gravity  is; 

w  ■  .961  w  . 
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